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[Article  by  N.  A.  Mukhstdlnov,  snginesr] 

[Text]  Rassarch  on  natural  convactlon  In  porous  aadla  with  a vertical 
teaverature  gradient  has  establlshad  that  air  (fluid)  aoveaent  occurs 
at  a cartain  critical  dlffersnce  in  teaveratures  [1,  2].  A definite 
velocity  and  reglaa  of  the  alr^ current  arise  in  a porous  aedium  for 
each  difference  in  teiqieracures  beyond  the  critical  point.  It  is  con- 
sidered that  the  laadnar  current  reglae  occurs  in  porous  nedla  if  the 
Reynolds  nuaiber  Re  <10,  while  the  regime  is  turbulent  if  Re  >7,000  [3]. 
In  the  range  of  Reynolds  nuaibers  Re  ” 10-7,000  air  aoveaent  corresponds 
to  what  is  callsd  the  transitional  current  regime  with  local,  not 
general,  turbulence  of  the  current  between  particles. 


He  know  that  for  the  laalnar  reglaa  of  an  air  current  (or  fluid)  losses 
of  pressure  are  proportional  to  velocity  in  the  first  degree,  while  for 
the  turbulent  reglae  it  is  the  second  degree. 


1.  Losses  of  Pressure  with  Air  Filtration  Through  Porous  Media  Under 
Conditions  of  Natural  Convection 


It  has  been  proven  both  theoretically  [4]  and  experiaantally  that  losses 
of  pressure  in  porous  media  for  any  Reynolds  nuiri»er  can  be  described  by 
the  following  two-tera  foraula: 


djp 

dy 


(1) 


where  dp/dy  are  pressure  losses,  is  the  modulus  of  d3maadc  viscosity 
of  the  air,  k and  b are  test  coefficients  which  depend  on  the  type  of 
porous  aadlua  and  the  physical  propertiss  of  the  air,  and  w is  the  rate 
of  air  filtration. 


In  foraula  (1)  the  first  tsra  characterises  the  influence  of  viscosity 
sad  the  second  is  the  influence  of  inertial  forces,  which  havs  bean 


studied  hy  many  investigators  [3,  5,  6].  A comparison  of  the  results  of 
different  investigations  shows  t^t  the  values  of  the  parameters  k and 
fluctuate  in  e broad  range.  At  the  saaw  time  to  receive  analytic  calcu- 
lation formulas  for  determining  pressure  losses  requires  the  introduction 
of  so  many  simplifying  assumptions  that  the  result  is  in  large  part  use- 
less [3].  In  calculations,  therefore,  empirical  and  seeiemplrlcal  expres- 
sions for  parameters  k and  b which  generalize  experimental  material  are 
usually  used.  An  acqualntanSe  with  existing  recooBendatlons  showed  that 
the  aforementioned  coefficients  have  not  yet  been  determined  for  condi- 
tions of  natural  convection  in  large-pore  media. 

We  carried  out  special  experimental  investigations  to  detenslne  pressure 
losses  and  the  filtration  regime  in  large-pore  media.  The  experiments 
were  conducted  under  laboratory  conditions  with  air  at  above-freezing 
temperatures. 

The  experimental  installation  was  a thermally  insulated  tube  (internal 
diameter  of  87  centimeters,  one  meter  long)  with  a steel  plate  bottom 
(five  Billllmeters  thick)  containing  20-mllllmeters  holes  evenly  distri- 
buted over  the  entire  area  of  the  bottom.  An  electric  heater  was  set  in 
the  lower  part  of  the  tube.  The  upper  part  of  the  unit  ended  in  an  exit 
cone  with  an  angle  of  90  degrees  at  the  apex  and  an  opening  9.5  centi- 
meters In  diameter.  For  the  experiment  the  tube  was  loaded  with  rocks 
and  thermocouples  were  set  in  the  pores  every  50  centimeters  (five  thermo- 
couples in  each  of  three  ranges),  then  the  exit  cone  was  secured  her- 
metically to  the  tube.  The  temperature  in  different  layers  of  the  tube 
was  recorded  by  the  thermocouples  with  a precision  of  0.3  degrees  C. 
while  air  velocity  at  the  exit  of  the  exit  cone  was  measured  by  a vane 
anemometer  with  a sensitivity  threshold  of  0.15  meters  per  second. 


Readings  were  taken  every  30  minutes.  The  air  In  the  base  of  the  tube 
was  heated  by  an  electric  heater.  Owing  to  the  uneven  distribution  of 
air  temperatures  in  the  porous  medium,  the  warm  air  rose,  giving  part  of 
its  heat  to  the  medium  loaded  In  the  tube  (the  rock) . 

The  electric  heater  switched  off  after  the  air  In  the  lowest  measurement 
range  reached  a temperatxire  of  100-120  degrees  C.  Continued  air  move- 
ment in  the  tube  occurred  entirely  through  difference  between  the  pore 
temperature  and  outside  temperature.  The  data  were  processed  on  the 
assumptions  that  air  was  evenly  distributed  in  the  tube  on  the  cross- 
section  and  that  the  lifting  force  was  proportional  to  pressure  losses 
at  any  moment  in  time.  The  letter  can  be  derived  from  the  following 
reasoning. 


Whan  sir  moves  under  the  influence  of  a lifting  force,  equation  (1)  for 
the  axis  colneidlag  with  the  direction  of  the  force  of  gravity  appears 
as  follows!  1.  4 

^w  + b^w^  -rPrviAT,  (2) 


where  T is  the  specific  weight  of  the  air,  ^ is  the  dilation  of  the 
elr,  a is  porosity,  and  ^ T is  the  temperature  difference  between 
pore  and  outside  eir. 


US  integrsts  equations  (1)  and  (2)  for  the  length  of  the  porous  layer 
In  the  experlaental  unit,  considering  parameters  s , k,  and  b_  and  the 
rate  of  filtration  w to  be  constant. 

As  a result  we  obtain  for  equation  (1) : 

(Pn  — (1 ' ) 

and  for  equation  (2) 

IPiu.  — Pmuj It. K -L -s)H  b.,w^H,  ^2 * ) 

where  is  the  air  pressure  at  the  input  of  the  tube,  pe^<  is  the 

outlet,  (Pet  n is  the  pressure  difference  with  forced 

convection,  i*  the  same  for  natural  (free)  convection,  and 

W is  the  elevation  of  the  layer  of  porous  medium. 

The  right  parts  of  equations  (1*)  and  (2')  can  be  made  equal  in  experi- 
ments by  selecting  the  corresponding  dlffsrenceof  pressures  or  lifting 
force.  In  this  case,  equating  the  right  and  left  parts  of  the  equations 
(1*)  and  (2*)  we  obtain: 


(Pn Pwuf)»t*u.  1C  — IPul  “jPowJc)c»  ■}* 


When  the  cause  of  air  movement  is  the  difference  in  temperatures  of  the 
pore  and  outside  air  the  condition  p^,  ^ p is  met  with  adequate  pre- 
cision. It  then  follows  from  equation  (3)  that  p p 

there  is  a formal  coincidence  of  the  pressure  gradient,  for  forced 
movement  with  a lifting  force  in  equation  (1)  and  for  natural  convection 
in  equation  (2). 


The  regime  of  an  air  (fluid)  current  in  a porous  medium  is  determined  by 
a diagram  of  the  relationship  between  the  coefficient  of  hydraulic  resis- 
tance in  filtration  (fp)  and  the  Reynolds  number  (7]. 


In  processing  the  results  of  our  tests  the  coefficient  of  resistance  (f_) 
and  the  Reynolds  number  (Re)  were  computed  according  to  formulas  [6,  7J: 


f-5/r.U  Ttfdm* 

“ (1  — ’ 

„ K'd 

“bav(l  - my 


(4) 

(5) 


***®>f*  t the  acceleration  of  the  force  of  gravity,  d is  the  effective 
dlmmeter  of  the  cleavages  of  the  porous  medium,  a.  i«  the  ratio  of  a 
cleavage  surface  to  the  surface  of  a sphere  with  diameter  d,  ^ is  the 
kinematic  viscosity  coefficient  of  the  air,  A T._  -0co-tHaD  the 
average  temperature  difference,  t„«p  is  the  outside  tehperafc  ure  of  the 

*****  *■*  ®***  temperature  in  the  pores  inside  the  exper- 

iSMntal  device  • 
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Flgur*  1.  Dlagri 


(1)  Haating  Coll; 

(2)  8Uu7: 

(3)  Hood  Chipa; 

(4)  Rocka; 

(5)  Thanocooplaa; 

(6)  Staal  Diaphraga  with  Holaa; 

(7)  Mirror  GalwanoMCar; 

(8)  Thoraootat. 

itaralaa  dc.  it  waa  aaauaad  that  tha  taavaratura  of  tha  air 
poraa  chuigaa  ■onltonleallyt  that  la*  saoothly.  In  tha 
;ura  waa  racordad  In  thraa  saetlona.  Baa ad  on  tha  thraa  nai 
)£  pora  taaparaturaa  a cnrra  (parabola)  waa  drawn  to  aatab- 
idanca  f(y)  ( 0 la  tha  air  tawparatura  In  tha  poraa  and 
llnata) . 6ct>  waa  dataxwlnad  by  tha  fonwU 


/(y)«/y 


whara  I la  tha  alawatlon  of  tha  layar  in  tha  anparlMntal  unit; 
a - 2(6i  ♦ 03>*4d2;  b - 482-01-303;  O3  la  awaraga  air  ta^»aratura  in  tha 


por«s  in  th«  lo««r  aeasureaent  section  of  the  tube;  62  Is  the  saiae  for 
the  Biddle  section;  6^^  Is  the  same  for  the  top  section. 

The  average  temperature  In  the  measurement  sections  was  determined  by 
readings  from  the  five  thermocouples.  The  readings  on  one  plane  dif- 
fered from  one  another » a result  of  the  uneven  movement  of  air  on  a cross- 
section  of  the  experimental  unit.  Losses  from  the  unit  at  the  Input  and 
outplet  were  not  taken  Into  account  In  procasslng  experimental  results 
(see  Figure  1).  One  of  the  tests  was  made  without  the  steel  diaphragm 
to  determine  the  effect  of  losses  at  the  input. 

The  diaphragm  was  replaced  by  a grid  with  cell  dioMnslons  of  20  x 20 
centimeters.  The  resiilts  showed  no  significant  differences.  In  all 
likelihood  losses  at  the  Input  do  not  exceed  the  precision  of  measure- 
ment of  losses  In  the  tests. 

Losses  at  the  outlet  were  estimated  by  calculation,  which  showed  that 
losses  at  the  outlet  were  less  than  the  lifting  force  contained  In  the 
exit  cone  2-3  hours  after  the  beginning  of  the  experiment.  The  tests 
were  made  with  stones  with  average  effective  diameters  of  10.6  and  19 
centimeters.  The  average  effective  diameters  were  determined  by  reduc- 
tion to  the  diameter  of  an  equal-sized  sphere.  The  porosity  of  the  medium 
In  the  tests  varied  from  0.42  to  0.47. 


The  tests  results  are  shown  In  Table  1 (partially)  and  In  Figure  2.  The 
ev*tage  values  of  test  results  In  the  range  of  Reynolds  number  from  10 
to  30  are  described  by  the  equation 


0.85. 


(7) 


In  Lev's  formula  [5]  for  determining  drop  of  pressure  In  a layer  the 
followltig  values  of  the  coefficient  of  hydraulic  resistance  are  recom- 
mended: for  Re  <10,  fp  - 100;  Re  - 50,  fp  - 5;  Re  - 100,  f - 1.2-1. 3; 

Re  > 1,000,  f^ 0.8-1. 2.  It  la  apparent  that  tha  value  ox  the  coeffici- 
ent of  resistance  for  Re  ■ 50  according  to  Lev's  data  Is  half  the  size 
of  the  same  figure  computed  according  to  formula  (7).  In  our  opinion, 
the  Increase  In  resistance  of  the  porous  medium  to  air  movement  under 
conditions  of  natural  convection  Is  caused  by  strong  cross-shifting  of 
flowing  air  becatiee  the  reason  for  the  moveswnt  Is  heat  exchange  between 
the  air  and  the  faces  of  the  cleavages.  At  some  moments  In  time  It  Is 
even  possible  for  a closed  circulation  to  form  in  certain  pores,  which 
reduces  the  live  sections  of  the  channels.  Moreoever,  errors  are  possi- 
ble in  the  teats  themselves.  First,  the  pores  which  we  selected  from 
the  whole  set  are  random  ones.  The  flow  of  air  In  these  pores  and  Its 
c<*P«ratttre  may  not  be  typical  of  all  phenomena  In  the  particular  section 
of  the  tube  at  any  particular  moment.  Second,  the  introduction  of  measur- 
ing Instruments  (anemometers,  mercury  thermometers)  Into  the  flow  disturbs 
the  established  regime  of  air  current  In  tha  tube.  Third,  the  porosity 
of  tha  madl«a  varied  In  the  tests  (from  0.42  to  0.47). 
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It  should  be  noted  that  using  ssiall 
values  for  the  coefficient  of  resis- 
tance In  calculating  the  thermal 
regime  of  the  lower  prism  of  rock- 
fill  dams  yields  rapid  changes  In 
the  temperature  of  the  prism  which 
are  not  observed  In  practice.  This 
Is  especially  true  when  the  coef- 
ficient of  turbulent  filtration 
(bg)  Is  small. 

The  coefficients  of  permeability 
(k)  and  turbulent  filtration  (bg) 
of  the  porous  medium  corresponding 
to  the  results  of  our  experiments 
are  determined  from  equation  (7) 
by  the  following  dependences: 


*^3100(1  -m)*’ 
b..  = 0.85 

Figure  2.  Dependence  of  Coefficient 
of  Resistance  of  Porous  Medium  on  These  are  correct  for  a porous  medium 
Reynolds  Number  (Test  Data)  effective  diameters  of  more  than 

10  centimeters.  Thus,  experimental 
Investigation  determined  the  coef- 
ficients needed  to  solve  the  problem  of  nonlinear  filtration  of  air  In 
large-pore  media. 


(1  -/«)?> 
dm?  ■ } 


2.  Statement  of  the  Problem  and  Method  of  Solving  It 

Work  [8]  obtains  an  equation  of  the  following  type  for  determining  air 
movement  In  porous  media  with  natural  convection: 


^ (k,  ^ (b 

( ’5x1  dy\  ' dy ) dx 


where  ^ Is  the  flow  function. 


k' 


(9) 


The  anss  velocities  along  the  coordinate  axes  are  expressed  through  the 
flow  function  by  the  formulas: 


(10) 


where  p is  air  density  and  u and  v are  a projection  of  the  velocity  vector 
on  axes  OX  and  OT. 


It  was  assumed  In  deriving  equation  (9)  that  pressure  losses  are  propor- 
tional to  mass  velocity  In  the  first  degree,  that  Is,  coefficient  k does 
not  depend  on  velocity.  It  can  be  shown  that  equation  (9)  Is  also  proper 
where  coefficient  fc^  Is  *a  function  of  velocity. 


In  the  latter  case,  in  conformity  with  [9]  the  coefficient  kj^  can  be 
determined,  with  due  regard  for  change  in  the  density  of  the  air,  from 
the  expression: 


*1  -/(pw) 


npii  Re  < Re.p-, 

Re>Re«i.. 


(11) 


where  Re  !•  the  critical  Reynolds  number  which  determines  the  passage 
of  the  air  current  from  one  regime  to  another. 

Considering  the  equation  for  the  average  temperature  of  the  cleavages 
of  the  fill  and  the  air  in  the  pores  we  obtain  the  following  system  of 
equations  for  determining  the  thermal  regime  of  the  lower  prism  of  a 
rock-fill  dam: 
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(T. 

^ (t  — 6V 

^ O',  fmdy  ox  rm  ox  oy  “ Co»m ' (13) 
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t): 


(12) 


where  t is  the  average  tenq>erature  of  the  cleavages  of  the  fill;  a is 
the  coefficient  of  temperature  conductivity  of  the  fill  (disregarding 
convection) ; is  the  volumetric  coefficient  of  heat  exchange;  Cg^  is 
the  volumetric  heat  capacity  of  the  fill;  c^J  is  the  volumetric  heat 
capacity  of  the  air;  t is  time. 

The  boundary  and  initial  conditions  of  equations  (12) , (13) , and  (14) 
were  discussed  in  detail  in  [8].  So  with  the  exception  of  certain  refine- 
ments we  will  not  cite  them  fully  here.  When  determining  the  heat  flow 
from  the  lower  prism  toward  the  line  of  its  contact  with  other  air- 
permeable  materials  a coefficient  of  transverse  convective  heat  conduct- 
ivity must  be  introduced  (by  analogy  with  forced  convection).  For  rocks 
in  natural  regimes  M.  E.  Aerov  and  N.  N.  Umnlk  consider  the  values  of 
the  coefficients  of  transverse  and  longitudinal  heat  conductivity  to  be 
close  and  propose  the  following  dependence  to  determine  it: 

X,  = + 5*>,  Re  Pr.  (15) 

where  App-  10.5  is  the  coefficient  of  convective  heat  conductivity  of 
the  air  with  i— in*r  current;  As  is  ths'^coefflclent  of  heat  conductivity 
of  the  air;  Pr  - i*  Prandtl  number;  Re  - 4w/F  is  the  Reynolds 

number;  F is  the  afea  of  the  surface  of  particles  in  a unit  of  volume  of 
the  fill;  B A.  is  the  experimental  constant. 

The  boundary  conditions  for  equations  (13)  and  (14)  (d^/ dn  ■ 0 and  8 *t|,ap) 
are  more  strictly  realized  at  4-5  meters  from  the  daylight  surface  of  the 


lower  prim  because  Che  Influence  of  Che  cemperacure  of  Che  air  moving 
ouC  of  Che  lower  prim  of  a rock>flll  dm  Is  llccle  felc  here. 


The  Cransfer  of  che  boundary  condiClons  of  equaClons  (13)  and  (14)  by 
Che  IndlcaCed  dlscance  corresponds  Co  InCroduclng*  on  che  permeable  layer 
of  che  lower  prim,  a flcClclous  porous  layer  wlch  permeablliCy  equal  Co 
chac  of  Che  fill  buc  wlch  zero  volumeCrlc  heac  capaclCy.  On  che  surface 
of  che  flcclclous  layer,  where  dd/dx  • 0,  Che  condlclon  df/dn  Is  scrlccly 
realized.  This  follows  from  Che  orchogonallcy  of  che  harmonic  funcclons 
which  saclsfy  equaclon  (14)  wlchouC  Che  rlghchand  pare. 


The  sysem  of  equaClons  (1)  was  used  Co  decermlne  che  chermal  regime  of 
a crlangularly  shaped  lower  prim  (Figure  3)  wlch  an  Inlclal  testperacure 
of  0 degrees  C.  The  boundary  condlclons  of  Che  problm  are  shown  In 
decall  in  Figure  3.  Ic  was  assumed  chac,  as  funcClons  of  coordlnaces, 
changes  In  Che  denslcy,  specific  welghc,  vlscoslcy,  and  dllaclon  of  che 
air  in  che  pores  would  noc  play  a signlflcanc  pare  In  decermlnlng  che 
velocicy.  _ 


Figure  3.  Olagrm  of  che  Calculaclon  of  che  Thermal 
Regime  of  a Triangular  Lower  Prim  (B  ■ 38m;  H <■  21  m) 


Then  equaClons  (13)  and  (14)  of  sysem  (1)  cake  che  form 

1 dC.  o6  •]  (7y  gfi  _r 

^ ‘ m ay  ox  m ox  ~ 


where  Po  vlscoslcy,  denslcy,  and  dllaclon  of  che 

air  corresponding  Co  che  average  ceoqteracure  of  che  air  moving  In  che 
pores. 


Th«  v«lu«  of  function  <^(z,  7t  '^ ) dlscrot*  point*  on  the  surface  of  Che 
fictitious  layer  was  detenined  froa  condition  d^/dn  ■ 0 according  to  the 
following  formula: 


where 


w Ax  sin  a.\.  Axsln'x,  //',  i.vsins..', 


The  solution  to  equation  (17)  was  considered  found  if  the  difference  of  the 
two  investigation  approziaationa  satisfied  the  condition: 


(24) 


where  € is  a given  comparison  constant. 


For  known  values  of  function  (^(z,  y* T ) at  discrete  points  In  the  domain 
under  consideration  it  is  necessary  to  determine  Che  flow  function  by  the 
following  equation: 


(•».  y.  -) 


(25) 


In  this  case  the  rate  of  filtration  w is  determined  by  the  formula 


Ue  approzlsMCe  the  derivatives  d^/dz  and  dp/dy  by  the  one-sided  spatial 
differences,  namely 


Wh  assuBM  that  ^Q.i,j  ^a,J-l  known  and  do  not  depend  on  the 
value  j.  Then,  substitutl^  (26)  into  (25)  we  perform  the  integration 
within  t&e  bound*  of  the  discrete  point  with  coordinates  z - n(Lz  and 
y • jAy  (1]]. 


As  a result  we  obtain  the  following  algebraic  equations: 
a)  for  c>0;A»  4ac-»)^>0 


brl  c V ri  / 4c 


1^0  (4iic  — Y-) 


. . Arsh 
bcl'c 


X,  . 


(27) 


b)  for  c > 0;  ^ • 0. 


- -?K,  + ^ <2»> 

wh«r«7  Is  tbs  kinsastlc  viscosity  coefficient  of  the  air  corresponding 
to  the^average  taaperature  of  the  air  Moving  In  the  pores 

iy»  )’ 

where  Arsh  Is  the  Inverse  value  of  the  main  hyperbolic  sine. 

To  slapllfy  the  calculation  of  foraulas  (27)  and  (28)  the  ratios  below 
■ay  be  used: 

Arsh  (—*)»—  Arsh  r, 

Arsh  (*)  - In  (r  - 1 r^i. 

In  this  case  ln(s  -f  ] s^  -f  1)  is  a natural  logarltha 

The  solution  to  equations  (27)  and  (28)  «ust  be  begttn  froa  the  lower  left 
comer.  In  this  case  ^q_i  4 and  ..2.  f'*  determined  froa  the  boundary 
conditions  and  do  not  In  fict  depenii^on  vq  j.  Thus,  the  values  of  0 *nd 
a when  datersdnlng  the  next  value  of  the  flow  function  at  a discrete 
point  are  always  definite.  The  solution  to  equations  (27)  and  (28)  can 
be  found  using  one  of  the  methods  of  approxlaate  solution  of  algebraic 
equations.  The  flow  functions  on  the  surface  of  the  fictitious  layer  are 
determined  by  equation  (23)  as  the  result  of  a corresponding  substitution 
of  (p(x,  y,T)  for  4>(x,  y,Tr). 

The  solution  for  linear  filtration  on  condition  b^  ■ o Is  deterained  from 
equations  (27)  and  (28)  directly. 


The  numerical  values  of  the  parameters  Included  in  the  calculation  depen- 
dences were  taken  to  be  the  following: 


d • 0.3  ■; 

■ - 0.35; 

fo  • 1/273  1/degree; 

B - 38  ■; 

A 2 

fc  - 4.5*10~®  m^; 
1 ^ b - o; 

bo/f>o  • 200  1/m; 
. ^offo  " !/•» 


II  /k  - 3-10‘®  m2; 

bo  - 0;  , 2, 

Y ■ 4.79*10"2  fVhour; 
g*  ■ 9.81  m/sec^; 

- 6.19‘10"^  kg/hr/a^ 


g*  ■ 9.81  m/sec^; 

- 6.19‘10"^  kg/hr/a^; 

C|,  ■ 0.24  gcal/kg 'degree; 

C«i  ■ 0.34  Bcal/a^* degree; 

Ag*  2.1'10~'  gcal/a'houT'degree; 
P - 17.5  m^; 


* Meters  per  hour  squared  were  used  for  calculating  acceleration  of  gravity. 


B\-  0.076; 
a ■ 0.002  a^/hour; 

Xcoap  ■ 0*7  gcal/a.hour‘dagraa; 
■ 360  gcal/a2.dagraa; 


t “ 0*  C.; 
t^p  - 40*  C.; 
? - 0.95; 

€ - 0.001. 


Tha  voluaatrlc  coafflcla&t  of  haac  axchanga  was  dataralnad  ac  dlscrate 
pointa  by  a calcu^tlon  according  to  the  foraula  _ 


whara  F is  tha  araa  of  tha  aurfaca  of  cleavagea  In  a unit  of  voluaa;  <p  la 
a coafflclant  which  takaa  account  of  the  decreaae  In  haat  axchanga  result- 
ing froB  tha  intamal  gradlant  of  taaq>eratures  in  tha  claawagasy 

is  tha  local  coafflclant  of  haat  axchanga;  Cp  is  tha  haat  capacity  of 
tha  air  at  constant  praasura;  w^  j Is  tha  rata  of  filtration  at  tha 
point  with  coordlnatas  x • n/kx  aM  y • fdy. 

Tha  tlaa  stap  (At)  during  calculations  was  variabla  and  was  determined 
froa  the  stability  of  tha  explicit  finite  difference  aquation  for  average 
teaperatura  of  tha  cleavages  of  tha  lower  prism  [8].  Tha  thermal  reglae 
of  tha  lower  priaa  was  calculated  in  two  variations  (first  variation  whan 
k “ 4.5*10~6  square  aetars»  1-bo  “ 0*  ^*'!{oPo  “ ^ aatar;  3-bgpn  “ 500 

1 aatar;  second  variation  whan  k • 3*10~‘^  square  aatars  and  bp  ■ 0).  All 
other  paraaatars  raaalnad  unchanged. 

Tha  results  of  tha  calculations  are  reprasantad  in  Figure  4 in  tha  fora 
of  graphs  of  tha  dapandanca  of  change  in  tha  average  taaparature  of  tha 
lower  prlsa  over  tlM.  With  a decrease  in  tha  paraeabillty  coefficient 
(k)  or  an  increase,  where  othor  conditions  are  equal,  in  tha  coefficient 
of  turbulent  filtration  (bp)  tha  rata  of  change  in  average  teaperatura  of 
tha  lower  prlsa  decreases.  Otherwise,  it  takas  twice  as  such  tiaa  to  cool 
tha  lower  priaa  to  tha  same  tac«paratura  where  bpp.  ■ 200  1 aatar  (k  ■ 
4.5*10~^  square  aatars)  ss  where  bp  “ 0.  This  is  related  to  tha  de- 
crease in  air  velocity  in  tha  lower  prlsa  of  tha  rock-fill  dam  (sea 
Figure  5).  Consideration  of  tha  nonlinearity  of  filtration  leads  to 
a aora  even  distribution  in  tips  of  tha  expenditure  of  air  passing  through 
tha  lower  priaa  (Figure  5).  It  should  be  noted  that  there  is  no  change  in 
the  qualitative  picture  of  change  in  the  teaperatura  of  the  lower  priaa 
for  any  coefficients  of  turbulent  filtration. 

Conclusions 

1.  Eapirlcal  dependences  were  obtained  for  deteralnlng  the  coefficients 
of  perasability  (k)  and  turbulent  filtration  (bp)  in  a large-pore  aediia. 


2.  A asthodology  is  proposed  for  nuaerical  solution  to  the  nonlinear  prob- 
lem of  air  filtration  in  an  isotropic  large-pore  aedlua. 

3.  Consideration  of  the  nonlinearity  of  air  filtration  in  pores  under 
conditions  of  natural  convection  decreases  the  rate  of  change  in  tha  teap- 
erature  of  the  lower  priaa  of  a rock-fill  daa. 


I 


Figure  4.  Chang*  In  Average  Teapereture  of  Lower  Prli 
(B  ■ 38  ■,  H ■ 21  b)  over  Tla* 


Figure  5.  Chang*  In  Filtration  Bat*  over  Tin*  at  the 
Creat  of  thn  Lower  Frias  at  Point  A 
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